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ABSTRACT 

A simple analytical model for estimating the fraction (f2g as ) of matter in gaseous form within the 
collapsed dark matter (DM) haloes is presented. The model is developed using (i) the Press-Schechter 
formalism to estimate the fraction of baryons in DM haloes, and (ii) the observational estimates of the 
star formation rate at different redshifts. The prediction for f2 gas from the model is in broad agreement 
with the observed abundance of the damped Lya systems. Furthermore, it can be used for estimating 
the circular velocities of the collapsed haloes at different redshifts, which could be compared with future 
observations. 

Subject headings: - cosmology: large-scale structure of universe - intergalactic medium - quasars: 
absorption lines 
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high column density systems are important in understand- 
ing the baryonic structure formation because they contain 
fair amou nt of the neutral hydrogen in the universe at high 



redshifts (Wolfe et al. 1986, Storric-Lombardi, McMahon 



& Irwin 1996). In contrast, the ionised hydrogen is mostly 



contained in the low column density systems, which mani- 
fest themselves as the Lya forest in the quasar absorpt ion 



spectrum (see, for example, Miralcla-Escudc et al. 1996 
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in the DLAs is similar to the baryonic mass which can 
collapse and efficiently cool in the DM haloes. However, 
not all the collapsed baryonic matter in the haloes remain 
in gaseous form - a fraction will turn into luminous stars. 
It is, therefore, clear that estimation of the total baryonic 
mass in the collapsed haloes and the mass turned into stars 
can be used to calculate the mass of the cooled gas in the 
haloes. 

In this paper, we use the simple idea that the mass con- 
tained in DLAs (i.e., the mass which remains in gaseous 
form in the collapsed haloes) should be equal to the to- 
tal baryonic mass in the collapsed haloes minus the mass 
which is turned into stars. In doing this, we have neglected 
the mass of the ionised gas within the haloes (which is 
small because the gas is dense enough to be shielded from 
the background radiation) and also the small amount of 
mass which might have turned into molecular gas. The 
total baryonic mass contained within the collapsed ha loes 
can be estimated through Press-Schechte r formalism ( Mc 
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Currently, the theoretical understanding of the DLAs 
i s based on the hierarchical models for ga laxy formation 
(|White fe Rees 1978| , |White fc Frenk 199l|) . These models 
start with the fact that, initially, the DM density inho- 
mogencities collapse via gravitational instability and form 
potential wells. In the next stage the baryonic matter, 
mostly in gaseous form, follow these DM potential wells. 
The baryonic gas can cool and form galaxies, provided 
its cooling timescale is short compared to the Hubble ex- 
pansion timescale. Once this condition is met, the gas is 
able to dissipate its energy and 'virialise' in the centres 
of the DM haloes. This cooled gas, contained within the 
DM haloes, is believed to be responsible for the damped 
Lya absorption lines in the quasar spectra. Various ana- 



lytical and semi analytical models (Mo fc Miralda-Escude 



ter which has turned into stars is calculated from the ob- 
served estim ates of the star form ation rate (SFR) at var- 
ious epochs ( Bteidel et al. 199£ ). It turns out that this 
simple model is able to produce the mass of the cooled gas 
in the collapsed haloes as seen in observations. Further, 
it is able to predict the circular velocities of the collapsed 
haloes associated with DLAs. 

This simple model should not be thought of as a com- 
petitor to the previous detailed models which use a lot 
more inputs for their modelling. Rather, we have tried to 
show that it is possible to reproduce many of the physi- 
cal effects through a simple argument of 'baryon conser- 
vation'. Our paper should be judged in this backdrop of 
simplicity. 

The paper is arranged as follows. The next section dis- 



1 



2 



cusses the basic analytical formalism used for the model. 
It describes the usage of the Press-Schechter formalism to 
calculate the total baryonic mass in collapsed haloes. It 
also discusses how to obtain the stellar mass density from 
the observed SFR. In section || we discuss the parameters 
of our model, followed by the results. In the last section, 
we summarise our main conclusions, and discuss the limi- 
tations of our model. 

2. ANALYTICAL FORMALISM 

This section contains the basic formalism for the model. 
Altho ugh moot of it io available in oxioting literature, 
briefly 



We will assume that all the baryons are either converted 
into luminous stars, or they form gaseous clouds. Mathe- 
matically, this can be expressed as 



n^(z) = fi Afiga .(z) + n,(z) 



(6) 



It is believed that the damped Lya systems are mainly 
contributed by the gaseous clouds within the collapsed and 
virialised dark matter haloes. Hence, one should not in- 
clude the low column density systems in this analysis, 
which are believed to be density p erturbations within a 



diffuse intcrgalactic medium (IGM) (Choudhury, Padman 
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the mass contained within the collapsed objects of a given 
mass range, we us e the Press Schechter formalism (Press 
& Schechter 1974 ). The total mass (dark matter) con- 



tained per unit comoving volume within collapsed objects 
of (logarithmic) mass range (InM, ln(M + dM)) is given 

by 



p$°{z) dlnM 
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Po- 



dlncr(M) 
dlnM 



(1) 



where po is the mean comoving density of the universe. 
The other quantities are defined as 



M = f 



3 (sin kR — kR cos kR) 



k 3 R 3 



fc 3 -P DM (fc) 
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The quantity v is defined as 
v{z,R) = 



D(z)a(R) 



(2) 



(3) 



ft)ivi(fc) is the normalised dark matter power spectrum, 
D(z) is the growth factor for density perturbations and S c 
is the critical density, usually 1.69 for VL m = 1 flat universe. 
The corresponding f2 is given by 
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Po 



(4) 



where p c = 2.8 x 10 n /i 2 M Mpc~ 3 is the present critical 
density of the universe and Q m — pa/p c - 

We shall now calculate the total mass contributed by 
baryons in the collapsed haloes within the logarithmic 
mass range. In order to do this, we assume that the bary- 
onic fraction of matter in each halo is same as the global 
value. Then the fl contributed by baryons in collapsed 
haloes within a logarithmic mass range (In M, ln(M+dA/)) 
is given by 
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^ M) ^AT nb (5) 



We now have the total baryonic ? contained within 
collapsed dark matter haloes of mass M at a given epoch. 
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It is clear that once we obtain the quantity f2» contained 
in stars, we can estimate QM.gas from the above relation. 
To calculate f2* from observational data, we proceed as 
follows: The star formation rate (SFR) (dp*/di) = p*(t) 
is defined as the rate at which baryonic mass is converted 
into stars per unit (comoving) volume. Given this quan- 
tity, we can obtain the total density of matter contained 
within stars at a particular t 



p*(t) = dt p\(t) 
Jo 



or, in terms of redshift 



p*(z) — [ dz p\(z)-^- 
The corresponding 17 is 



P*(z) 



H(z)(l + z) 



£l*(z) 
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One can determine p*(t) from observations, which can 
then be integrated to give f2» (z). It is well known that not 
all the gas which is turned into stars is removed from the 
gaseous phase forever. Actually, some of the stellar ma- 
terial is returned by stellar winds and supernovae. How- 
ever, we shall ignore this contribution for the following 
reasons: (i) It is difficult to implement this feedback effect 
i n our model wit hout introducing more free parameters 
( Efstathiou 2000| ). This spoils the simplicity as well as 
the predictive capacity of our model, (ii) The SFR has 
large uncertainties at high redshifts due to extinction (see 
below). The correction due the feedback through stellar 
winds and supernovae is less than this uncertainty in the 
SFR, and can be ignored as a first approximation. 

To compare our results with observations, or to discuss 
any observational consequences of our model, it is better 
to work in terms of the circular velocities (v c ) of the col- 
lapsed haloes, rather than in terms of their masses. The 
mass M and v c can be related to each other using spherical 
collapse model. The releva nt equations for a universe w ith 
cosmological constant are (Bomerville fc Primack 1999) 



gm n A m rf„ 



and 
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2GM 



(10) 
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where A v i r (z) is the virial density of the collapsed halo at 
redshift z and r vlr is the corresponding virial radius. The 
circular velocity is calculated assuming that the virialised 
halo has a singular isothermal density profile p(r) oc r~ 2 . 
The above relations can be used to eliminate r v ; r and ob- 
tain 



of v c . It is clear that haloes of lower circular velocities 
contain more mass at higher redshifts. 

There is quite a bit of uncertainty about the SFR. This 
is mainly because of our lack of knowledge of the extinction 
due to dust. The extinct ion has been modelled in various 
ways (see, for example, Stcidcl ct aL 19991, lAdelbergei 
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et al. 2001). In this work, we take the SF R data ob- 



H 2 (z)A M .{z) 



taine d from ultraviolet luminosity observations ( Lilly et al 



aH-2( y )A vir ( y ) 



Thus, for a given z we can relate M to v c , provided we 
know A vir (z). This result can then be used for obtaining 
the mass density, either dark matter (0^ al ° (2;)) or baryonic 
i^v^b ( z ))> contributed by haloes having a circular velocity 
v c . 

Given the background cosmological model, our formal- 
ism can now be used for estimating the quantity Q Va ,g&s(z) 
for different values of v c , which can then be compared with 
observations. Note that we have only one free parameter 
in our model, namely v c . It turns out that we can, in 
principle, predict the circular velocities of collapsed haloes 
at different redshifts by comparing our model with DLA 
observations. 

3. MODEL PARAMETERS AND RESULTS 

We have considered two different cosmological models 
with the parameters listed below: 

SCDM: n m = 1, A = 0, h = 0.65, n b h 2 = 0.02 
LCDM: n m = 0.3, Q A = 0.7, h = 0.65, Vt b h 2 = 0.02 
The critical density for collapse S c is only weakly depen- 
dent on the ft a for flat cosmological models (Eke, Cole, & 
Frcnk 1996| ); hence we take it to be 1.69 for both the mod- 
els! The next cosmological input that is required is the 
form of the DM power spectrum. We take the fol lowing 
form for Pdm(&) flEfstathiou, Bond, fc White 1992[) 



Ak 



(1 + [ah + (/3k) 1 - 5 + (7k) 2 ]") 2 /- 



(13) 



where v = 1.13, a = {SA/T)^ 1 Mpc, (3 = 
{S.O/T)^ 1 Mpc, 7 = (1.7 '/^h' 1 Mpc and T = Q m h. The 
normalisation parameter A is fixed through the value of 
(is = o~( R = 8/i -1 Mpc) We take the values of as to be 
given by (|Eke et al. 1996; ) 



0"8 



(0.52 ± o.o4)ji-°- 46+olon ™ (if r> A = 0) 
(0.52 ± 0.04)si-°- 52+013Om (if n A = 1 - n m ) 

(14) 

Given the above parameters, we can calculate the Press- 
Schechter mass function. However, since we prefer to work 
in terms of v c rather than M, we need to know the quan- 
tity A V i r (;z) (see equation (Il2)). For flat models with a 
cosmological constant, t his is given by the fitting formula 
([Bryan fc Norman 1998|) 



A vir (z) = 187r 2 + 82a;-39a; 2 ; x = Q m (z) - 1 (15) 

In Fig. |l|, we plot the quantity rijj al ^(z) for the two dif- 
ferent cosmological models and for three different values 



1996|, |Connolly et al. 1997] , |Madau et al. 1996|, gteide 



et al. 199S), corrected for extinction by (steidel et al 



1999). The data can be fitted with a function of the form 



a e 



(16) 



e cz + d 

with the parameters given by 

a = 0.13 M e yr- 1 Mpc" 3 , b = 2.2, c = 2.2, d = 6.0 (17) 

There are two effects which may modify the form and 
the absolute value of this SFR. The first modification is 
due to cosmology. The above data is for a flat il m = 1 
universe, with h = 0.5. In order to use the same points 
for a different cosmology, one has to note that the SFR 
is proportional to luminosity per volume, which in turn is 
inversely proportional to the distance for a given redshift. 
This means, while considering a flat universe with a cos- 
mological constant (O m + Qa = 1) we have to multiply the 
above SFR by a correction factor 



Ccosm (z) 



h 

05 



2[l-(l + z)- 1 /2] 



f z dz< [(l-O m )+ft m (l + z') 3 ]- 1/2 



(18) 

In Fig U we plot this correction factor for the two cosmo- 
logical models of our interest. The correction for SCDM 
model comes just because we are using h = 0.65, and is 
independent of z. 

The second effect is due to the stellar initial mass func- 
tion (IMF). The SFR is usually calculated from the ob- 
served luminosity using some particular IMF. The SFR 
data we are using is calculated using the standard Salpeter 
IMF wit h a slope T = —1.3 5 for a mass range I OOMq 



O.LMq ( [Madau et al. 1996| , |Stcidel ct al. 1999[ ). How- 
ever, the real IM F probably becomes flatter below 1M@. 
Lcithcrcr (1998| ) calculates that the extrapolation of the 
Salpeter IMF to O.IM© overestimates the SFR by a factor 
of 2.5. In this work, we take into account the flattening 
and eventual turning over of the IMF below IMq by sim- 
ply reducing the SFR by a factor 2.5. 

Implementing both the corrections, we get the corrected 
SFR in the form 



a e bz C c , 



2.5 



(19) 



with the parameters a, b, c, d given in equation (|17|). Given 
P*(z), we can calculate the quantity £l*(z) using equations 
(||) and (||). This quantity is plotted in Fig. |^, for two cos- 
mological models. 

We now obtain the main result of our work, ^l Vc ,gas(z), 
and compare it with DLA observations. We have al- 
ready discussed how to compute the quantities £i^ al £(z) 
(see equation (|j)) and SI* (see equations (||) and ([)])). 
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Fig. 2. — The correction factor for the observed SFR because of using a different cosmological model. The observed data points are for a 
flat Q m = 1, h = 0.5 universe, while we are interested in cosmological models with different parameters (see text for the cosmological models 
considered in this paper). The correction factor is plotted as a function of z for the two cosmological models of our interest (SCDM and 
LCDM). The solid line denotes the correction for SCDM model (arising because we use h = 0.65), while the dashed one is for the LCDM 
model. 
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Fig. 4. — The density of matter contained in gaseous form n„ c> g a s(z) as a function of z, plotted for two cosmological models. The circular 
velocity ranges from v c = 200 km s _1 to v c = 250 km s _1 in both the plots. 
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Given these two quantities, we can analytically compute 
r2 1 , e gas (z) using equation (||). The observed data used 
for comparing our model is obtained from 



Peroux et al 



(20010. First we consider a scenario in which the circular 
velocities of the collapsed haloes do not have significant 
evolution, i.e., they are more or less constant over the red- 
shift range. To be more precise, we consider a velocity 
range from v c = 200 km s _1 to v c = 250 km s _1 for both 
the cosmological models and see how they compare with 
the observed points. Fig. || shows the comparison within 
this velocity range. As can be seen from the figure, the 
ball park estimate of Q Vc ,gas(z) from our simple model is 
well within observational constraints. Furthermore, the 



generi d trend of evolution of Q Vc ^ as (z) is also reproduced. 

In passing, we mention that we have plotted the mass 
density r2„ cjgas (z) contributed by haloes within a logarith- 
mic velocity range. We could have as well plotted the 
quantity Q gas (z), which is the mass density contributed 
by haloes having a wide range of circular velocities ( say. 



from v r 



50 km s to v c 



250 km s" 1 ) (see 



& Mil alda-Escude (1994) )) . Mathematically, this quantity 
can be obtained from fl Vt ,^ s (z) through the relation 



^gas (-0 



-ft 



v c ,gas 



d In v r 



-m- 



However, when compared with observations, this quan- 
tity turns out to be much higher than the observed mass 
density contributed by DLAs. This means that either (i) 
most of the observed DLAs are hosted by DM haloes with 
a small range of circular velocities or, (ii) we are severely 
underestimating the SFR. At this stage, it is quite difficult 
to distinguish between the two possibilities. 

Let us now return back to Fig. ||. One notices that ob- 
servational points in higher redshifts are better fitted by 
curves with lower velocities. For example, for the LCDM 
model, the point at z = 2.35 (fourth from left) has a bet- 
ter match with the curve having v c = 250 km s _1 (a curve 
with a higher v c will do still better), while the point at the 
highest redshift [z — 3.92) can be fitted with a curve of 
lower velocity (v c < 200 km s _1 ). This might motivate one 
to consider a second scenario where the circular velocity 
falls as redshift increases. 

It should be obvious that it is impossible to constrain 
the evolutionary pattern of v c with such large error bars in 
the observational data. However, we shall illustrate how it 
can be done in principle. We consider a simple evolution 
of v c parametrised as 



Vc(z) 



250 km s- 



= V 



1 



(21) 



where 77 is a parameter of the order unity, to be fixed by 
comparing with observations. Fig. [5] shows the compari- 
son between our model and observations when v c is given 
by the above relation. The values of r\ for the two cosmo- 
logical models are 



V 



1.167 (for SCDM) 
1.163 (for LCDM) 



(22) 



The simple evolution law of v c matches with observa- 
tions reasonably well, especially at z > 2. The match is 



not that good for z < 1, and one has to introduce more 
complicated redshift dependence to take care of this. We 
do not perform such exercise in this work, as it is difficult 
to constrain any parameter with such large observational 
errors. 

We are thus able to show that our model, in spite of 
being extremely simple, is able to reproduce the broad ob- 
servational trends for the quantity 0„ cigas (z) (the density 
of baryons contained within collapsed haloes in gaseous 
form). The model can, in principle, be used to predict the 
circular velocities of the collapsed haloes as a function of 
z. 

It should be noted that the detailed model of Kauffmanr] 
(19961 ) gives somewhat lower circular velocities compared 
to ours. At this stage it is difficult to co mpare our results 



directly with those of [Kauffmann (1996 ). Since both have 
used the PS to calculate the mass in collapsed objects, 
the difference must lie in the estimation of SFR. We have 
used an observed estimate for the SFR. It is not obvious 



from the work of Kauffmann 1996 how the results for SFR 



match with observations. 

As an extension of this work, we also calculate the dis- 
tribution of circular velocit ies of the DM haloes, using the 
Press - Schechter formalism flHaehnelt, Steinmctz, fc Rauch 



1998, Hachnclt, Steinmctz, k, Rauch 2000). However, we 



preferred to keep this paper simple and focused on one key 
idea only. 

4. DISCUSSION 

We have presented a simple analytical model for esti- 
mating amount of matter left in gaseous form within the 
collapsed DM haloes. The only ingredients for this model 
are (i) the Press-Schechter formalism for the collapse of 
DM haloes, and (ii) the observational estimates of the SFR 
at different redshifts. Our model reproduces the mass con- 
tained in the baryonic gas which has cooled within the 
DM haloes as expected in observations. Furthermore, the 
broad trends seen in the observed data are reproduced in 
our model. We have used it to estimate the circular ve- 
locities of the collapsed haloes at different redshifts, which 
might be compared with future observations. 

The model has a couple of important uncertainties, de- 
scribed below: 

(i) The observational estimate of the SFR is quite uncer- 
tain, mainly because of the extinction. There are different 
models for the extinction in the literature, and the result- 



ing estimate diff er quite a bit (Steidel et al. 199E, Hop- 



kins et al. 2001). It is difficult to prove what is the best 
method to use at high redshift for estimating the extinc- 
tion, since there are so few constraints. In our model we 



have used the extinction corrected data of Steidel et al, 
|(1999| ), as th eir estimates arc consistent with the sub- mm 
background (Adelberger & Steidel 2000) 



Note that this uncertainty may have a large influence 
on the values of v c cal culated in this paper. For example, 
the SFR estimated by Hopkins et al. (2001) is about 5 to 
6 times higher than what we have used at z > 2.5. Using 
such a SFR would have lead us to use a higher value of 
baryonic matter density so as to match the observations 
of Sl Vcygas (z). The only obvious way to increase the total 
baryonic matter is to include the haloes with lower ve- 
locities. Thus there is a possibility that the haloes with 
lower v c than what is calculated in this paper might also 
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Fig. 5. — The density of matter contained in gaseous form Vc:gllB (z) as a function of z, plotted for two cosmological models. The circular 
velocity is assumed to be proportional to the scale factor. It is evident from the figure that this z-dependence of v c does not hold for z < 1. 



contri bute to f^ eigas (z). Thus, the absolute values of v c 



predicted in this work might have errors. However, the 
broad features (like v c decreases with increasing z) are ex- 
pected to be still valid. 

(ii) In our model, we have neglected the feedback from 
supernovae and stellar winds because there is no simple 
way to deal with them. This might have resulted in a 
slight underestimate of ri 1 , cjgas (z). However, we do not 
expect this correction to modify the broad trends of the 
model. 

We have found that our model matches the observations 
in the redshift range 2 < z < 5, provided the circular ve- 
locities of the halo es are within the range 100 km s -1 < ■ 
v c < 250 km s" 1 . |Haehnelt, Steinmctz, fc Rauch (1998|) 



and Hachnelt, Steinmctz, & Rauch (2000) co mpare hy 



drodynamical simulations and observed data (Prochaska 



& Wolfe 1997) to conclude that the typical circular veloc- 
ities of the haloes at high redshifts are ~ 100-200 km s _1 . 
This value is in agreement with our analysis. Also, kine- 
matic studies of Lcdoux ct al. (1998) find that there is a 
trend for the velocity to decrease with redshift. The circu- 
lar velocity cannot be predicted for lower redshifts (z < 2) 
because of the large errors on the observed Q Vc ,gns(z). 
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